Decoherence processes in a current biased dc SQUID 
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A current bias dc SQUID behaves as an anharmonic quantum oscillator controlled by a bias current 
and an applied magnetic flux. We consider here its two level limit consisting of the two lower energy 
states |0) and |1). We have measured energy relaxation times and microwave absorption for different 
bias currents and fluxes in the low microwave power limit. Decoherence times are extracted. The 
low frequency flux and current noise have been measured independently by analyzing the probability 
of current switching from the superconducting to the finite voltage state, as a function of applied 
flux. The high frequency part of the current noise is derived from the electromagnetic environment 
of the circuit. The decoherence of this quantum circuit can be fully accounted by these current and 
flux noise sources. 

PACS numbers: Valid PACS appear here 



In the past years, coherent manipulation of two and 
multi-level quantum systems, efficient quantum read- 
outs, entanglement between quantum bits have been 
achieved [3 13 El demonstrating the full potential 
of quantum logic in solid state physics. At present, future 
developments require longer coherence times 0, 0, ITof . 
In contrast with atomic system, the huge number of 
degree of freedom makes its optimization a challenging 
problem. Up to now, the most successful strategy has 
been to manipulate the quantum system at particular 
working points where its coupling to external noise is 
minimal 2]. Nevertheless, experimental analysis of deco- 
herence phenomena in superconducting circuits remains 
a priority for a full-control of quantum experiments. Dif- 
ferent models for the noise sources have been proposed 
to describe the decoherence processes acting on various 
qubits0|lS[l3. However a complete and consistent un- 
derstanding of decoherence remains a current and open 
problem. In this paper, we study decoherence processes 
of a phase qubit: the current biased dc SQUID. 

This superconducting circuit consists of two Joseph- 
son junctions (JJ), each with a critical current Iq and a 
capacitance Cq. The junctions are embedded in a super- 
conducting loop of inductance L s , threaded by a flux <Iv 
In the limit where L s Io w $o/2ir, the phase dynamics of 
the two junctions can be mapped onto a fictitious particle 
following a one dimensional path in a 2D-potential Q. If 
the biasing current If, is smaller than the SQUID criti- 
cal current I c , the particle is trapped in a cubic poten- 
tial well characterized by its bottom frequency u> p (Ib, $&) 
and a barrier height AU(Ib, $&) (Fig.n] a )- The quantum 
states in this anharmonic potential are denoted \n), with 
corresponding energies E n , n — 0, f , ... In the following, 
only the lowest states |0) and |1) will be involved. For 7f, 
well below I c , these two levels are stable and constitute 
a phase qubit. 

When the bias current is close to 7 C , AC/ de- 
creases and becomes of the order of a few hu; p . The 
ground state can tunnel through the potential barrier 
and the SQUID switches to a voltage state[ll|. The 
tunnelling rate Tq of the ground state |0) is given by 



the well-known MQT formula for underdamped Jj|l2fl: 
ro(ib, <£■;,) = aujp exp(— 36AU/5fkj p ), where a is of order 
unity. 

The environment of the dc SQUID induces fluctuations 
of the bias current and the bias flux. In this work, we 
show how the current and flux noise sources can be sep- 
arately quantified. This is achieved by escape measure- 
ments of the SQUID at specific working points where it 
is mostly sensitive to current or flux noise. Using these 
identified noise sources, the measured decoherence times 
are fit precisely as a function of bias current. 

Experimental results are analyzed assuming a linear 
coupling between the SQUID and the environment de- 
grees of freedom. We suppose that current SI and 
flux 5$ noises are generated by independent gaussian 
sources. Here, Sx (x = I or is an operator acting 
on the environment. Their fluctuations are specified by 
the quantum spectral densities 5vfi/)[l4j. In presence 
of flux microwave (MW) excitation, the total Hamilto- 
nian H in the SQUID eigenstates basis {|0) , |1)} reads: 
H = —^hvoiaz — Hvr cos(2irvt)a x + N where a x and a z 
are Pauli matrices and v m = (E\ — E )/h. The first term 
is the qubit Hamiltonian and the second term describes 
the MW excitation of reduced amplitude Vr at frequency 
v. In this notation, i/r is also the Rabi precession fre- 
quency for a tuned excitation (y — vqi) . The last term 
is the coupling to the noise sources. For our circuit it is, 
within linear response, 
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where 77 is the asymmetry inductance parameter (see 
below), rj(0) = cos 8 + 77 sin 0, r$(9) = sin 8 and 8 
is the angle between the escape and the mean slope 
directions in the 2D potential [ill ITo| . To first order, 
the transverse noise proportional to a x only induces 
depolarization. The longitudinal term proportional to 
a z induces "pure" dephasing. The qubit sensitivity to 
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FIG. 1: (a) Squid cubic-quadratic potential, (b) Electri- 
cal environment of the SQUID, (c) Calculated effective real 
impedance i? e // versus frequency. 



longitudinal noise is given by the partial derivatives 
(dvoi /dlb) and (dfoi/<9$&). They depend strongly on 
the experimental working point and increase near the 
critical current. 

The measured SQUID consists of two large aluminum 
JJs of 15 /im 2 area (J = 1.242 and C = 0.56 pF) 
enclosing a 350/xm 2 -area superconducting loop. The two 
SQUID branches of inductances L\ and contribute to 
the total loop inductance L s = 280 pH with the asymme- 
try parameter r\ = (Li — L2) / L s = 0.414. The immediate 
electromagnetic environment of the SQUID is designed to 
decouple the circuit from the external world. It consists 
of two cascaded LC filters (see Fig. \T\h). A large on- 
chip inductance L oc = 9 nH is made of two long and thin 
superconducting wires which value, derived from the nor- 
mal state resistance, is dominated by the kinetic induc- 
tance. The gold thin film parallel capacitor, C g ss 150pF, 
introduces a finite resistor. Its dc value at 30 mK is 
R g = 0.1 Vt giving the gold resistivity p g = 1.210 -8 fim. 
The second filter consists of the bounding wires, with an 
estimated inductance Lf = 3nH, and a surface mounted 
(SMC) capacitor C SM c = 2nF and four 500 SMC resis- 
tors. The nominal room temperature microwave signal 
is guided by 50 f2 coaxial lines, attenuated at low tem- 
perature before reaching the SQUID through a mutual 
inductance M s = 1.3 pH. Special care was taken in mag- 
netic shielding and bias lines filtering. All these electrical 
parameters were determined independently |lll ll3T | . Our 
environment model predicts two resonances (resistance 
dips in Fig. ^c). They were observed in a similar set-up 
and the associated resonance frequencies were in precise 
agreement with the model. 

The current noise through the SQUID comes mostly 
from its immediate environment thermalized at T = 
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FIG. 2: (a) Measured escape current (dots) versus external 
applied flux fitted to MQT theory (solid line) at 30 mK. (b) 
The width of the probability distribution P eS c(h) (dots) fitted 
to the 2D MQT predictions. The solid curve takes the low 
frequency flux noise into account while the dashed line does 
not. At bias flux $i,o (resp. $i,i) the sensitivity to flux noise 
is zero (resp. small) while it is maximum at $52 • 



30 mK {y T = k B T/h = 600 MHz < v QX ). The quantum 
spectral density of the current noise, Sj{v) in this en- 
vironment is set by the fluctuation-dissipation theorem: 

Sl(v) = Mcoth(^) + lJ-ReffM" 1 where ^cff^r 1 

is the real part of the environment circuit admittance. 
R e s(v) is calculated using the electrical circuit shown in 
Fig. Eh and is plotted in Fig. He To a good approxima- 
tion, the root mean square (RMS) current fluctuations 
are of order yj ksT/ L oc = 6 nA. Most of the noise is 
peaked around 30 MHz, a frequency much smaller than 
Vt- A simple estimate of the flux noise produced by the 
inductive coupling to the 50 £1 coaxial line shows it can 
be neglected in the following. 

The escape probability P CS c{h) out of the supercon- 
ducting states is measured at fixed flux using dc current 
pulses with At = 50 /xs duration and h amplitude. Each 
measurement involves 5000 identical current pulses and 
the total acquisition time is T m = 10 s. The escape cur- 
rent /esc is defined as the current lb where the escape 
probability P eS c{Ib) = 0.5 and the width of the switching 
curve A J = Ih — h as the difference between the currents 
where P esc (I h ) = 0.9 and P esc {h) = 0.1. In Fig. the 
dependence of I csc and AI on $b are plotted. By fit- 
ting the escape current curve I csc {^b), the experimental 
parameters of the SQUID (Iq, Co, L s , rf) are determined. 

Moreover, escape measurements are a sensitive tool to 
characterize noise (frequency range and amplitude). If 
noise frequencies exceed the inverse of a current pulse 
duration At^ 1 , the tunnel rate fluctuates during each 
current pulse. The escape probability is controlled by 
the average (Fo) escape rate in the frequency window 

^M~LM- p -c = 1 - exp [ - (r (/ 6 + SI, $ b + 6$)) At] 
nil hg. The current noise produced by the electrical 
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environment lies in this frequency interval. Its effect is 
to decrease I eac ($b) by about 6 nA, the RMS current 
fluctuations (unobservable in Fig. 0a). Similarly, the 
width of the switching curve is not affected. 

On the other hand, if noise frequencies are slower than 
At -1 , the tunnel rate is constant during a pulse, but 
fluctuates from pulse to pulse. In this limit, the escape 
probability becomes P csc = ( 1 — cxp [ — T (/& + SI, $fc + 
<5<i>)Ai]), where the statistical average () is now in fre- 
quency range from T~ to At . To first order, low fre- 
quency noise does not affect I csc , but increases the width 
AI. Thus AI is the best quantity to probe the origin 
and the magnitude of the low frequency fluctuations: if 
the flux <I>f, is set at the value <&bo which maximizes I c , 
the SQUID is only sensitive to current fluctuations since 
= 0. In the vicinity of this flux, the measured width 
is explained by the usual MQT theory. Hence the mea- 
sured RMS current fluctuations in the [T^ 1 , Ai _1 ] in- 
terval (low frequency current noise) is below 0.5 nA, the 
error bar in AI measurements. This is consistent with 
the 0.1 nA RMS value derived from the spectral density 
of noise at frequencies below At -1 . For other applied 
fluxes, the width is slightly larger than MQT prediction, 
indicating a residual low frequency flux noise. The de- 
pendence of AI on $b shown in Fig. 0b is perfectly ex- 
plained by a gaussian low frequency flux noise. Its RMS 

1 /2 

amplitude, (S&j^p) ~ 5.5 x 10~ 4< J>o, is extracted from 
the fit shown in Fig. [21b and is attributed to the flux 
noise in the [100 mHz, 20 kHz] frequency interval. The 
origin of flux noise may be due to vortices trapped in the 
four aluminum contact pads located at a 0.5mm distance 
from the SQUID. 

Hereafter we discuss dephasing and relaxation induced 
by the noise sources previously identified. These incoher- 
ent processes are experimentaly studied with low power 
spectroscopy and energy relaxation measurement. As de- 
scribed in Rcf . 5] , a MW flux pulse is applied followed by 
a 2 ns duration dc flux pulse to perform a fast but adia- 
batic measurement of the quantum state of the SQUID 
(Fig. 0c inset). The duration T M w = 300 ns of MW 
pulses is sufficient to reach the stationnary state where 
the population pi of the level |1) only depends on v and 
the amplitude vr. The microwave amplitude vr is cal- 
ibrated using Rabi like oscillations 0- In the two level 
experiments discussed in this paper, the measured es- 
cape probability P csc induced by the dc flux pulse can 
be interpreted as P OS c = fJsc + (Pell - P&H) x Pi{v, v R ). 
Pise denotes the escape probability out of the pure state 
\n). In Fig. |31 a and 0b, the escape probability versus 
microwave frequency v are plotted at two different bias- 
ing points. The experimental curves present a resonant 
peak which position and full width at half maximum de- 
fine the resonant frequency vq\ and Av. Spectroscopy 
experiments are performed in the linear regime and Av 
is experimentally checked to be independent of the MW 
amplitude. Relaxation measurements were performed by 
populating the |1) state with low power MW tuned at vq\ 
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FIG. 3: (a) and (b) Escape probability versus applied mi- 
crowave frequency with amplitude vn < 5 MHz at two dif- 
ferent working points (I b = 2.288 ^iA, $ b i = —0.117 <&o) and 
(h = 0.946 AtA,*M = -0.368 $ ), respectively. The points 
are experimental data and the continuous lines are the Fourier 
transforms of f co h(t) (see text), (c) Measured escape probabil- 
ity versus delay time (dots) fitted to an exponential law with 
Ti = 95 ns (continuous line). The inset specifies the timing 
of the measurement pulse which follows the MW excitation 
pulse. 



during a time Tmw = 300 ns, and measuring its popula- 
tion with increasing time delay Tdoiay after the end of the 
MW pulse. As shown in Fig. 0c, the escape probability 
follows an exponential relaxation with a characteristic 
time T\. In Fig. measured resonant frequency vqi, 
relaxation time Ti and the inverse of microwave width 
Av^ 1 are plotted versus current bias for the two differ- 
ent applied fluxes $bi (close to ${,o) an d $t 2 shown in 
Fig. Vqi, Ti and Au^ 1 decreases as J& gets closer to 
I c . For these two applied fluxes, the vm dependence fits 
perfectly the semiclassical formulas for a cubic potential 
|17j using the same SQUID electrical parameters as those 
extracted from escape measurements. 

The depolarization rate Tj -1 is given by the sum T-j -1 = 
Ffl+r^ of the relaxation Tr and the excitation T e rates. 
These two rates are calculated using Fermi golden rule. 
At low temperature, excitation can be neglected and Tr 
reads: 



r, 



4C hvo: 
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Neglecting the high frequency part of the flux noise, 
one obtains T\ = 2i? e ff (ym)CQ j r\{6) where R c fi{voi) — 
(2-kL oc Vqi) 2 / R g (vQi). R g (v) — aR s is the high frequency 
resistance of the gold capacitor where R s = ^Wp^pgD is 
the surface resistance and a is a dimensionless geometri- 
cal parameter. The T\ versus /{, dependence is well fitted 
with a = 200 as the only adjustable parameter (Fig. 0b). 
This value is the right order of magnitude for the geom- 
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FIG. 4: Resonant transition frequency (a), relaxation time (b) 
and microwave width (c) as function of bias current at $ti = 
—0.117 $o and $(,2 = —0.368 <E>o respectively right and left 
curves. Symbols correspond to experiments and continuous 
line to model predictions. 

etry of the gold capacitor. 

Relaxation alone is too weak to explain the value of Av 
and "pure" dephasing also contribute to the linewidths. 
First, we consider the time evolution of the reduced den- 
sity matrix in the basis {|0) , |1)} in the absence of MW. 
The linear coupling to noise sources induces a time de- 
cay of the amplitude f co h{t) of the coherence terms. Since 
current and flux noises are independent, f co h(t) is factor- 
ized as f co h(t) = fi(t) /*(*) exp(-2t/T x ) where //(£) and 
/$(t) are respectively the "pure" dephasing contributions 
due to current and flux noises. 

The current contribution is given by the well-known 
gaussian noise formula 

EEl: //(*) = exp[- \t 2 x 



(2tt^) 2 /+~ di/S/(i/)sinc(7n/i)], where {dv m /dl b ) is 
extracted directly from the slope of the experimental 
curve of Fig. 0Ja. We neglect flux noise contributions 
with frequencies higher than 20kHz. Since the acquisition 
time of absorption spectra and escape measurements are 
similar, the SQUID undergoes the same RMS flux fluctu- 
ations in the two experiments. In these conditions, 
takes the simple gaussian form: /$(t) = exp [ — ^t 2 x 

(2tt^) 2 ((5$| /F )] . (dvoi/d$ b ) was computed using the 
known electrical parameters of the SQUID |l9j. 

Within linear response, the shape of the resonance 
curve is proportional to the Fourier transform (FT) of 
fcoh(t). Resonance curves in Fig. [3]a and [3Jb are fit- 
ted using P esc - P\sl oc FT{f coh }(v - i/ 01 ) (continuous 
line). Our model explains perfectly the shape of the ex- 
perimental curves. In Fig.0]c, the theoretical width Av 
extracted from the curve FT{/ co ^}(i/) is in very good 
agreement with experimental data without free param- 
eter. When If, gets close to I c , the partial derivatives 
{dvQi/dlb) and (dvoi/d$b) increase: the noise sensitivity 
increases and Av broadens. For bias points correspond- 
ing to $(,2, the width is due to current and flux noise. 
For a bias flux equal to $bi, the effect of flux noise is 
small and the width is dominated by current noise. At 
this flux, for < 1.95 fxA, our model predicts satellite 
resonances around vqi which are not observed. Other 
noise mechanism may blur the predicted satellite peaks. 

In conclusion, we have shown how the flux and cur- 
rent noise present in this controlled quantum circuit can 
be separately identified. We measured the decoherence 
times at low microwave power where the quantum cir- 
cuit can be reduced to a two level system. Analyzing the 
coupling of the SQUID to the known noise sources, the 
measured relaxation times and the resonance width can 
be fully understood. 
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